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Abstract 


ShipMo3D  is  an  object-oriented  library  with  associated  user  applications  for  predict¬ 
ing  ship  motions  in  calm  water  and  in  waves.  This  report  describes  improvements 
to  the  treatment  of  hull  maneuvering  forces  and  rudder  forces  within  ShipMo3D.  A 
proportional- integral-derivative  (PID)  autopilot  is  also  introduced  to  ShipMo3D.  Pre¬ 
dicted  turning  circles  for  the  tanker  Esso  Osaka  show  good  agreement  with  results 
from  sea  trials.  Maneuvering  predictions  based  on  experimental  hull  maneuvering 
force  coefficients  are  better  than  predictions  based  on  hull  coefficients  from  a  regres¬ 
sion  method. 


Resume 


Le  ShipMo3D  est  une  bibliotheque  de  programmes-objets  munie  d’applications  utili- 
sateur  connexes  qui  permettent  de  predire  les  mouvements  d’un  navire  en  eaux  calmes 
et  dans  les  vagues.  Le  present  rapport  decrit  les  ameliorations  effectuees  an  traitement 
des  forces  de  manoeuvre  de  la  coque  et  des  forces  du  gouvernail  du  ShipMo3D.  Le 
ShipMo3D  est  egalement  muni  d’un  autopilote  proportionnel,  integral,  derive  (PID). 
Les  cercles  de  giration  predits  pour  le  petrolier  Esso  Osaka  demontrent  de  bonnes 
similarites  avec  les  resultats  des  essais  en  mer.  Les  predictions  qui  reposent  sur  les 
coefficients  de  force  de  manoeuvre  des  coques  experimcntales  sont  mcilleures  que  les 
predictions  qui  s’appuient  sur  les  coefficients  de  coque  de  la  methode  de  regression. 
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Executive  summary 

Improved  Maneuvering  Forces  and  Autopilot  Modelling 
for  the  ShipMo3D  Ship  Motion  Library 

Kevin  McTaggart;  DRDC  Atlantic  TM  2008-162;  Defence  R&D  Canada  -  Atlantic; 
September  2008. 

Introduction:  ShipMo3D  is  an  object-oriented  library  with  associated  user  applica¬ 
tions  for  predicting  ship  motions  in  calm  water  and  in  waves.  Motion  predictions  are 
available  in  both  the  frequency  domain  and  the  time  domain.  For  predictions  in  the 
frequency  domain,  a  ship  is  assumed  to  travel  with  quasi-steady  speed  and  heading  in 
waves.  For  predictions  in  the  time  domain,  the  ship  can  be  freely  maneuvering  in  ei¬ 
ther  calm  water  or  in  waves.  This  report  describes  improvements  to  evaluation  of  hull 
maneuvering  forces  and  rudder  forces.  The  report  also  describes  the  implementation 
of  a  proportional-integral-derivative  (PID)  autopilot  within  ShipMo3D. 

Principal  Results:  ShipMo3D  now  has  an  improved  method  for  evaluation  of  hull 
maneuvering  forces  based  on  measured  hull  maneuvering  force  coefficients.  Also,  pre¬ 
dicted  rudder  forces  now  duely  consider  incident  flow  velocities  and  resulting  drag 
and  lift  force  components.  Predicted  turning  circles  for  the  tanker  Esso  Osaka  show 
good  agreement  with  results  from  sea  trials.  Maneuvering  predictions  based  on  ex¬ 
perimental  hull  maneuvering  force  coefficients  are  better  than  predictions  based  on 
hull  coefficients  from  a  regression  method. 

Significance  of  Results:  This  report  concludes  the  development  of  Version  1.0 
of  ShipMo3D  and  addresses  the  final  concerns  that  arose  during  development.  It  is 
expected  that  future  revisions  to  Version  1.0  of  ShipMo3D  will  be  limited  to  minor 
bug  Exes. 

Future  Plans:  Development  of  a  new  version  of  ShipMo3D  is  progressing,  and  in¬ 
cludes  modelling  of  azimuthing  propellers. 
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Improved  Maneuvering  Forces  and  Autopilot  Modelling 
for  the  ShipMo3D  Ship  Motion  Library 

Kevin  McTaggart ;  DRDC  Atlantic  TM  2008-162 ;  R  &  D  pour  la  defense  Canada  - 
Atlantique ;  septembre  2008. 

Introduction  :  Le  ShipMo3D  est  une  bibliotheque  de  programmes-objets  munie 
d’applications  utilisateur  connexes  qui  permettent  de  predire  les  mouvements  d’un 
navire  en  eaux  calmes  et  dans  les  vagues.  Les  predictions  de  mouvement  sont  dis- 
poniblcs  dans  le  domaine  des  frequences  et  le  domaine  temporel.  Dans  lc  cas  des 
predictions  dans  le  domaine  des  frequences,  il  est  entendu  que  le  navire  navigue  a  une 
vitesse  quasi  constante  et  qu’il  met  le  cap  sur  les  vagues.  Dans  le  cas  des  predictions 
dans  lc  domaine  temporel,  le  navire  peut  manoeuvrer  librement  en  eaux  calmes  ou 
dans  les  vagues.  Le  present  rapport  decrit  les  ameliorations  apportees  a  L evaluation 
des  forces  de  manoeuvre  de  la  coque  et  des  forces  du  gouvernail.  Le  rapport  decrit 
egalement  Lajout  d’un  autopilote  proportionnel,  integral  et  derive  (PID)  au  sein  du 
ShipMo3D. 

Resultats  principaux  :  Le  ShipMo3D  est  maintenant  muni  d’une  methode  devaluation 
amelioree  pour  les  forces  de  manoeuvre  de  la  coque  qui  s’appuie  sur  des  coefficients 
de  force  mesures.  Par  ailleurs,  les  predictions  des  forces  du  gouvernail  prennent  main- 
tenant  en  consideration  les  vitesses  de  courant  incident,  la  resistance  obtenue  et  les 
composants  de  la  force  de  portance.  Les  predictions  des  cercles  de  giration  du  petrolier 
Esso  Osaka  demontrent  de  bonnes  similarites  avec  les  resultats  obtenus  lors  des  es- 
sais  en  mer.  Les  predictions  qui  reposent  sur  les  coefficients  de  force  de  manoeuvre 
des  coques  experimentales  sont  mcilleures  que  les  predictions  qui  s’appuient  sur  les 
coefficients  de  coque  de  la  methode  de  regression. 

Importance  des  resultats  :  Le  present  rapport  constitue  Letape  finale  de  l’elaboration 
de  la  version  1.0  du  ShipMo3D  et  aborde  les  dernieres  questions  qui  se  sont  presentees 
au  cours  de  celle-ci.  II  est  a  prevoir  que  les  prochaines  revisions  de  la  version  1.0  du 
ShipMo3D  se  limiteront  a  de  simples  corrections  de  bogues. 

Travaux  ulterieurs  prevus  :  L’elaboration  d’une  nouvcllc  version  du  ShipMo3D 
est  en  cours  et  comprend  la  modelisation  des  helices  en  azimut. 
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1  Introduction 


ShipMo3D  is  an  object-oriented  library  with  user  applications  for  predicting  ship 
motions  in  calm  water  and  in  waves.  Initial  ShipMo3D  developments  [1,  2,  3]  were 
aimed  at  predicting  motions  in  waves  for  a  ship  travelling  with  quasi-steady  speed 
and  heading.  ShipMo3D  was  subsequently  extended  to  predict  motions  in  calm  water 
and  in  waves  for  a  freely  maneuvering  ship  [4] .  This  report  presents  improvements  to 
maneuvering  forces  introduced  in  Reference  4.  A  revised  treatment  is  given  for  utiliza¬ 
tion  of  hull  maneuvering  coefficients  from  experiments.  Also,  a  revised  treatment  of 
rudder  forces  is  given,  which  includes  changes  to  rudder-propeller  interaction  forces. 
The  report  also  introduces  a  proportional-integral-derivative  (PID)  autopilot  model. 
The  improvements  described  in  this  report  are  included  in  ShipMo3D  1.0  [5,  6,  7]. 

This  report  includes  both  relevant  theory  and  validation  results.  Section  2  gives  axes 
systems  and  equations  of  motion  used  by  ShipMo3D.  Sections  3  and  4  present  im¬ 
proved  treatments  of  hull  forces  and  rudder  forces.  The  PID  autopilot  is  presented  in 
Section  5.  Section  6  gives  validation  of  turning  circle  predictions  for  the  tanker  Esso 
Osaka,  and  is  followed  by  conclusions  in  Section  7. 

2  Axis  Systems  and  Equations  of  Motion  for 
ShipMo3D  Computations 


Various  axis  systems  are  used  when  evaluating  the  forces  and  resulting  motions  for 
a  freely  maneuvering  ship.  Ultimately,  the  motions  of  a  freely  maneuvering  ship  are 
given  in  an  earth-fixed  axis  system  shown  in  Figure  1.  The  instantaneous  heading 
(to)  of  the  ship  is  denoted  x  and  the  heading  (from)  of  ocean  waves  is  denoted  v. 

Like  many  seakeeping  programs,  ShipMo3D  evaluates  ship  forces  and  resulting  mo¬ 
tions  in  translating  earth  axes,  as  shown  in  Figure  2.  When  ocean  waves  are  present, 
their  direction  relative  to  the  ship  is  given  according  to  Figure  3.  For  large  ampli¬ 
tude  motions,  forces  on  the  ship  can  be  evaluated  using  the  ship-based  axis  system 
of  Figure  4  and  then  transformed  to  translating  earth  axes. 

When  computing  ship  motions  in  the  time  domain,  it  is  necessary  to  select  a  suitable 
nominal  heading  and  nominal  speed  for  orienting  the  translating  earth  coordinate  sys¬ 
tem  used  to  solve  the  ship  motions  at  each  time  step.  Computations  with  ShipMo3D 
indicate  that  the  translating  earth  coordinate  system  can  be  oriented  based  on  the 
instantaneous  heading  and  ship  speed  along  that  heading  at  each  time  step. 

During  each  time  step  of  a  time  domain  simulation,  the  forces  and  resulting  accelera¬ 
tions  are  computed  in  translating  earth  axes.  The  resulting  accelerations  in  earth-fixed 
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Ship 


Figure  1:  Earth-Fixed  Coordinate  System 


Figure  2:  Translating  Earth  Coordinate  System 
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axes  are: 


(1) 

(2) 

(3) 


x*  =  f) i  cos x  —  Vr2  sin x 
Vf  =  ~Vi  sin  x  +  V2X 
X  =  -Ve 

Heave,  roll,  and  pitch  are  equivalent  in  earth-fixed  and  translating  axis  systems. 


The  equations  of  motion  in  the  time  domain  are  very  similar  to  those  presented  in 
Reference  4.  Changes  in  equations  of  motion  in  the  current  report  are  due  to  an  inte¬ 
grated  approach  being  used  for  evaluation  of  rudder  force  components.  This  integrated 
approach  for  rudder  force  components  is  required  to  obtain  accurate  predictions  of 
the  significant  forces  that  can  act  on  a  rudder  in  the  presence  of  a  propeller.  The  new 
equations  of  motion  are: 


m\  +  m  <*>)])  m)}  = 

-  B(U, oo,ij)l  {!)(«)}  -  [  [KhM(U,t-T)\  W(t)}  <It 

-  'c(U,oo^)]{V(t)}  +  {F\t)}  +  {FD(t)}  +  {Fcross(t)} 
+  {Fresist(t)}  +  {Fprop{t)}  +  {Frudder(t)} 


(4) 


The  ship  inertia  matrix  [M]  is  based  on  the  dry  ship.  The  added  mass  matrix 
[A{U,  ex:)]  is  the  infinite  frequency  value  dependent  on  ship  speed  U,  and  includes 
contributions  from  the  hull  and  appendages  as  follows: 

N f  oil  Nrudder 

[A(U,  oo)]  =  [Ahull(U,  oo)]  +  J]  [Af^]  +  [Arudder~i] 

i=l  i= 1 


Nskeg  Nbk 

+  [ A Ske9~i]  + 

i=l  i=  1 


(5) 


where  [Ahul1]  is  the  hull  added  mass,  is  the  added  mass  for  static  foil  i  (e.g., 

propeller  shaft  bracket),  ^Arudder~l]  is  the  added  mass  for  rudder  i,  [ Aske9~l ]  is  the 
added  mass  for  skeg  i,  and  [Abk~l]  is  the  added  mass  for  bilge  keel  i.  The  appendage 
added  masses  are  independent  of  speed  and  frequency  of  motion.  The  damping  matrix 
includes  contributions  from  the  hull  and  all  appendages  expect  for  the  rudders  as 
follows: 


B(U,oo,ri) 

N foil 

+  [Bhull-visc(U:\f)4\)\  +  [Bfoil-\U,  |t)4|)] 

i=l 

Nskeg 

+  Y  f/,1^1)]  +  ^[Bk-‘((7,|^4|)]  (6) 

i= 1  2—1 


_  '^hull—rad^jj-y  _|_ 


B 


hull— man 


(u,  it) 
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where  ^phull~rad]  is  the  hull  frequency-independent  damping  due  to  potential  flow, 
^hull-man ]  -g  (qie  hull  maneuvering  damping,  which  can  include  forces  from  lift  and 
cross-flow  drag,  [. Bhull~msc ]  is  the  hull  viscous  damping  due  to  roll,  [Bdml~l\  is  the 
damping  for  static  foil  i,  \^Bske9~l]  is  the  damping  for  skeg  i,  and  [Bbk~l]  is  the  damping 
for  bilge  keel  i.  The  retardation  function  matrix  \Khul1]  is  due  to  the  dependence 
of  hull  hydrodynamic  forces  on  frequency  of  motion.  The  stiffness  matrix  includes 
contributions  from  buoyancy,  hull  potential  flow  forces,  and  hull  and  appendage  lift 
forces  as  follows: 

C(U,  oo ,  ?7 )  1  =  [Cbuoy]  +  +  \chull~man{U,fi) 

N foil  N skeg 

+  E 

i=  1  i= 1 

Nbk 

+  £  [C7“_i(C0]  (7) 

i=  1 

where  [Cbuoy]  is  the  hull  buoyancy  stiffness,  \chul1  is  the  hull  frequency  independent 
stiffness  due  to  potential  flow,  \fjhull~man]  js  the  hull  stiffness  due  to  maneuvering 
forces,  [Cf°d~l]  is  the  lift  stiffness  for  static  foil  i,  [ Cske9~ ®]  is  the  lift  stiffness  for  skeg 
i,  and  [Cbk~l]  is  the  lift  stiffness  for  bilge  keel  i.  The  incident  wave  excitation  vector 
{F1}  consists  of  the  following: 

Nf0ii 

{F^t)}  =  {FI~hull(t)}  +  {F1-10^^)} 

1=1 

Nskeg  Ngk 

+  (8) 

i= 1  i=  1 

where  {Fr-huU}  is  the  incident  wave  force  on  the  hull,  is  the  incident  wave 

force  on  foil  i,  js  the  incident  wave  force  on  skeg  i,  and  is  the 

incident  wave  force  on  bilge  keel  i.  The  wave  diffraction  force  vector  {FD}  typically 
only  consists  of  hull  forces,  unless  panelled  appendages  are  included  in  diffraction 
computations  [3]. 

Reference  4  describes  computation  of  terms  {i?cross|)  j presist  {ppropy  arising  from 
hull  cross-flow  drag,  resistance,  and  propulsion.  The  term  {Frudder}  represents  all 
forces  acting  on  the  rudders,  with  the  exception  of  added  mass  forces  that  are  eval¬ 
uated  on  the  left  hand  side  of  Equation  (4).  Forces  in  { Frudder }  arise  from  lift  and 
drag,  and  include  effects  from  rudder  deflection,  nearby  propellers,  and  water  motion 
in  the  seaway.  Section  4  describes  the  evaluation  of  these  forces. 

Some  ship  force  components  are  most  easily  computed  using  ship-fixed  axes,  followed 
by  translation  to  translating  earth  axes.  Figure  4  shows  ship-based  axes,  which  can 


-Cfo*-i(U)\  +  ^  [Cske9~\Uj 
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vary  significantly  from  translating  earth  axes  when  motion  amplitudes  are  large. 
Within  ShipMo3D,  ship  resistance  and  propulsion  forces  are  assumed  to  act  along  the 
xs  direction  of  ship-based  axes.  If  buoyancy  and  incident  wave  forces  are  computed 
using  integration  of  panel  pressures  in  the  time  domain,  then  these  components  are 
also  evaluated  using  ship-fixed  axes. 

After  selected  force  components  are  computed  in  ship-based  axes,  they  are  trans¬ 
formed  to  translating  earth  axes  using  the  following: 


Fi(t) 

=  F\(t)  COS  7/5  COS  7/6 

—  F^it)  (-  cos  7/4  sin  t/6  +  sin  7/4  sin  775  cos  t)q) 

+  F^(t)  (cos  7/4  sin  7/5  COS  7/6  +  sin  7/4  sin  7/6 ) 

(9) 

m) 

=  F\  (t)  COS  7/5  sin  7/6 

+  F^  (t)  (cos  7/4  COS  7/6  +  Sm  7/4  sin  7/5  sin  7/6 ) 

+  Fg  (t)  (cos  7/4  sin  7/5  sin  7/6  —  sin  /74  COS  T)q) 

(10) 

W) 

=  —F±(t)  sin  7/5  +  F.^t)  sin  7/4  cos  7/5 

+  F^(t)  cos  7/4  cos  7/5 

(11) 

W) 

=  Ff(t) 

(12) 

m) 

=  Fg(t)  cos  7/4  -  Fq  (t)  sin  7/4 

(13) 

m) 

=  F^(t)  sin  7/4  +  F$  (t)  cos  7/4 

(14) 

where  Fj  is  the  force  for  mode  j  in  translating  earth  axes,  and  F^  is  the  force  for 
mode  j  in  ship-based  axes. 

For  a  ship  executing  a  turning  circle,  several  terms  are  used  to  describe  the  turning 
circle  characteristics,  as  illustrated  in  Figure  5.  Dimensions  are  typically  referenced 
to  the  point  at  which  rudder  deflection  was  initiated.  The  advance  refers  to  the 
displacement  in  the  direction  of  the  initial  heading,  and  the  transfer  refers  to  the 
displacement  lateral  to  the  the  direction  of  the  initial  heading.  The  terms  advance 
and  transfer  are  both  often  applied  to  either  any  point  on  a  turning  circle  or  to  the 
specific  point  at  which  the  ship  heading  has  changed  by  90  degrees  relative  to  the 
original  heading.  The  tactical  diameter  is  the  transfer  when  the  heading  has  changed 
by  180  degrees. 

3  Evaluation  of  ShipMo3D  Hull  Maneuvering 
Forces  from  Experimental  Coefficients 


This  section  provides  an  update  to  hull  maneuvering  coefficients  introduced  in  Refer¬ 
ence  4.  ShipMo3D  hull  maneuvering  coefficients  do  not  include  forces  due  to  potential 
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flow;  thus,  ShipMo3D  subtracts  potential  flow  forces  from  experimentally  obtained 
hull  maneuvering  coefficients.  In  Reference  4,  it  was  assumed  that  experimental  hull 
maneuvering  coefficients  are  measured  in  translating  earth  axes  at  high  motion  fre¬ 
quencies.  In  reality,  experimental  hull  maneuvering  coefficients  are  typically  measured 
in  stability  axes  at  low  motion  frequencies.  Note  that  stability  axes  follow  the  yaw 
motions  of  the  ship,  as  discussed  in  Schmitke  [8]. 

The  formulation  for  maneuvering  coefficients  presented  previously  in  Reference  4 
assumed  that  the  ship  drift  angle  was  small,  and  the  ship  forward  speed  was  approx¬ 
imately  the  same  as  the  total  ship  velocity  in  the  horizontal  plane.  It  is  no  longer 
assumed  that  drift  velocity  is  small,  and  ShipMo3D  hull  maneuvering  coefficients  are 
now  based  on  total  ship  velocity  in  the  horizontal  plane,  given  by: 


Note  that  total  ship  velocity  is  commonly  used  when  defining  hull  force  coefficients  in 
the  maneuvering  literature  [9,  10].  ShipMo3D  hull  maneuvering  coefficients  for  sway 
and  yaw  are  consequently  now  expressed  as  follows: 


td  hull— man 
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The  non-dimensional  sway  velocity  v'  is  given  by: 


The  non-dimensional  yaw  rate  is  given  by: 


he  L 
V 


(20) 

(21) 


ShipMo3D  hull  maneuvering  force  coefficients  are  evaluated  using  non-dimensional 
force  coefficients  based  on  Inoue  et  al.  [10].  The  linear  force  coefficients  in  Equations 
16  to  19  are  given  by: 


TDhull—man  ^  „  r  rp  \/ /  ahull— rad— stab  t  ,  n)  /oo\ 

-°22V  =  ~2PL  lmid  Yy  ~  B 22U  e  =  OJ  (22  j 

TDhull—man  ^  „  r  2  rp  \r /  Tjhull— rad— stab /  ,  n't  /no) 

&26V  =  ~2PL  lmidYr  ~  B2QU  \Ue  =  Uj  (26) 

TDhull—man  ^  „  r  2  rp  at/  ahull— rad—  stab  <  ,  n)  z'O/1) 

^62V  =  ~2PL  1  mUl  Nv  ~  =  Uj  (24) 

TDhull—man  ^  „  r  3  rp  AT1  id  hull— rad/ ,  ,  rA  /'oci'i 

Bmv  =  ~2pL  1  mid  7V?'  -  ^ee u  (u}e  =  0)  (25 j 
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where  p  is  water  density,  L  is  ship  length  between  perpendiculars,  and  Tmid  is  draft  at 
midships.  The  non-dimensional  force  coefficients  Y/n  Y/,  N'v,  and  N'r  include  both  vis¬ 
cous  and  potential  flow  effects.  The  terms  B^~rad~stab ,  B^l}~rad~stab ,  B^~ra d~stab, 
and  B^^~rad~stab  are  potential  flow  terms  in  stability  axes,  which  are  included  in  ex¬ 
perimentally  measured  force  coefficients.  These  terms  are  subtracted  when  obtaining 
maneuvering  force  coefficients  for  ShipMo3D  because  ShipMo3D  already  considers 
the  potential  flow  forces.  The  hull  radiation  damping  terms  to  be  subtracted  from 
the  above  equations  are  based  on  coefficients  presented  in  Reference  2,  with  terms 
here  given  for  stability  axes: 

' — "  hull  , 

dAjk  (U  =  0,  cae  =  0) 
dx 

derivative  of  potentials  giving  added  mass 
and  ue  =  0  denoting  zero  encounter  fre¬ 
quency. 


ryhull— rad—  stab f ,  ,  n\ 

ujkU  K  =  uj  = 

- — -  hull 

where  dAjk  /dx  is  determined  from  the  x 
with  U  =  0  denoting  zero  forward  speed 


4  Rudder  Forces 


This  section  gives  revisions  to  the  treatment  of  rudder  forces  presented  previously 
in  References  3  and  4.  These  revisions  are  based  partly  on  material  presented  in 
recent  books  by  Lewandowski  [11]  and  Faltinsen  [12],  Several  aspects  of  modelling 
of  rudder  forces  have  been  revised.  ShipMo3D  previously  considered  rudder  forces  to 
have  separate  force  contributions  from  ship  forward  speed  and  the  flow  induced  by 
the  propeller.  This  previous  treatment  has  been  replaced  by  an  integrated  approach 
that  considers  both  the  ship  forward  speed  and  the  flow  induced  by  the  propeller 
when  evaluating  the  flow  incident  to  the  rudder  and  the  resulting  total  force.  When 
evaluating  the  lateral  flow  incident  to  the  rudder,  ShipMo3D  can  now  include  the 
flow  straightening  effects  due  to  the  presence  of  the  hull  and  propellers,  which  can 
cause  a  significant  reduction  in  the  lateral  velocity  of  incident  flow.  Another  significant 
revision  is  the  treatment  of  rudder  lift  and  drag  forces.  It  was  previously  assumed  that 
the  total  force  acting  on  the  rudder  was  due  to  lift  forces  acting  normal  to  the  rudder. 
This  approach  is  used  in  other  manuevering  models,  such  as  those  in  References  9  and 
13.  ShipMo3D  now  uses  the  approach  described  by  Newman  [14],  whereby  lift  forces 
act  perpendicular  to  incident  rudder  flow  and  drag  forces  act  parallel  to  incident 
rudder  flow.  Limits  on  lift  and  drag  forces  at  high  angles  of  attack  are  also  now 
included  in  ShipMo3D. 

4.1  Lift  and  Drag  Forces  on  Rudders 

Lift  and  drag  forces  acting  on  a  rudder  are  shown  in  Figure  6.  The  lift  force  acts 
perpendicular  to  the  incident  flow  Vr,  while  the  drag  force  acts  in  the  same  direction 
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as  the  incident  flow  Vr.  The  lift  force  is  given  by: 

=  \psv?  c“'‘(a)  (27) 

where  S  is  the  rudder  area,  Vr  is  the  incident  flow  velocity,  Ch-’t  is  the  lift  force 
coefficient,  and  a  is  angle  of  attack  between  the  incident  flow  and  the  rudder.  At 
small  angles  of  attack,  the  lift  curve  slope  dCh^t /da  is  linear  and  can  be  evaluated 
based  on  Whicker  and  Fehlncr  [15],  as  discussed  in  Reference  3.  To  account  for  the 
influence  of  stall  effects,  ShipMo3D  places  a  limit  on  the  maximum  lift  coefficient  of 
a  rudder.  Data  presented  in  Reference  16  indicate  that  a  value  of  1.2  is  a  suitable 
nominal  limit  for  the  absolute  value  of  the  lift  coefficient. 


rprudder 

“lift 


Figure  6:  Lift  and  Drag  Forces  Acting  on  a  Deflected  Rudder  in  Incident  Flow 
The  drag  force  acting  on  the  rudder  is  given  by: 

F£it  =  \pSV?C*«>  (28) 


Bertram  [16]  and  Lewandowski  [11]  indicate  that  the  drag  coefficient  is  proportional 
to  a2  at  low  angles  of  attack,  and  can  be  represented  as  follows: 


(jdrag 


Qfjdrag 

d(a 2) 


2  dCdra 9 
d(a 2) 

'dclift\2  i 

da  J  nae 


(29) 

(30) 
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For  a  rudder  truncated  at  the  hull,  the  effective  aspect  ratio  ae  is  equal  to  2 s/c,  where 
s  is  the  rudder  span  and  c  is  the  mean  chord  length.  To  model  the  influence  of  large 
angles  of  attack,  ShipMo3D  applies  a  limiting  value  to  the  drag  coefficient  of  a  rudder. 
A  suggested  limiting  value  for  the  drag  coefficient  can  be  based  on  flow  normal  to  a 
flat  plate,  which  has  a  drag  coefficient  of  1.17. 


4.2  Rudder  Incident  Flow  Velocities  and  Resulting 
Forces 

The  most  challenging  aspect  of  computing  rudder  forces  is  the  evaluation  of  the 
incident  flow  into  the  rudder.  When  evaluating  the  incident  flow,  it  is  essential  to 
consider  ship  motions,  the  presence  of  the  ship  hull,  and  the  presence  of  propellers. 

For  the  present  analysis,  rudder  velocity  components  ur ,  vr,  and  wr  will  be  consid¬ 
ered,  which  are  velocities  of  the  rudder  relative  to  the  surrounding  fluid.  The  signs  of 
these  components  are  the  same  as  for  the  x,  y,  and  z  coordinates  in  the  translating 
earth  axes  of  Figure  2.  For  simplicity,  the  present  discussion  assumes  small  ship  mo¬ 
tion  angular  displacements;  however,  computations  within  ShipMo3D  evaluate  rudder 
forces  in  ship-based  axes  and  transform  them  to  translating  earth  axes.  Unlike  much 
of  the  maneuvering  literature,  the  present  discussion  considers  rudders  with  arbitrary 
orientation  rather  than  just  conventional  rudders  that  are  oriented  downward  from 
the  hull. 


The  effective  longitudinal  inflow  velocity  considers  the  presence  of  propellers,  and  is 
evaluated  using  the  following  equation  based  on  Faltinsen  [12]: 


where  Nprop  is  the  number  of  propellers,  ~pr op  is  the  rudder-propeller  interaction 

coefficient  from  propeller  j,  U  is  the  ship  forward  speed,  wr  is  the  rudder  wake  fraction, 
and  u siip—j  is  the  velocity  of  the  propeller  slipstream  at  the  rudder.  The  rudder- 
propeller  coefficient  Crudder~prop  will  have  a  value  approximately  equal  to  the  fraction 
of  rudder  area  within  the  propeller  slipstream.  The  velocity  of  the  propeller  slipstream 
can  be  evaluated  from  Bertram  [16]  based  on  original  work  by  Soding  [17,  18].  The 
flow  velocity  far  aft  of  the  propeller  based  on  potential  flow  is  given  by: 

U<slip  —  U  (1  — ■  Wprop )  y/l  +  Cth  (32) 

where  wprop  is  the  propeller  wake  fraction.  The  propeller  thrust  loading  coefficient 
Cth  is  given  by: 


jpprop 

1/2  p  U2  (1  -  Wprop)2  vr/4  D2 
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where  Dprop  is  the  propeller  diameter.  The  theoretical  slipstream  radius  far  aft  of  the 
propeller  is: 


r 


OO 

slip 


'  prop 


\ 


(34) 


where  rprop  is  the  propeller  radius.  The  slipstream  radius  based  on  potential  flow 
varies  with  longitudinal  location  as  follows: 


pot 
^ slip 
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0.14  (r siip /r prop)  +  r<slip/ ^ prop  X  [(^; 
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Similarly  the  slipstream  velocity  based  on  potential  flow  varies  with  longitudinal 
location  as  follows: 


u 


pot 

slip 


(x) 


Uslip 


(36) 


Due  to  turbulence  effects,  corrections  to  both  the  slipstream  radius  and  velocity  can 
be  made  as  follows: 


T*  slipi^C) 


^ slip( pt) 


pot 
^ slip 
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<l(x)  -  U  ( 1 
n%{x)  +  U  (1 
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—  w, 


prop ) 


Wprop  ) 

^  +  U  (1  —  Wprop )  (37) 


The  lateral  and  vertical  velocity  components  of  the  rudder  are  evaluated  as  follows: 


Vr  =  f]2  +  Xu  ft  f]6  (38) 

wr  =  /)3  -  xi% ft  775  (39) 


When  evaluating  the  effective  incident  flow  acting  on  a  rudder,  a  flow  straightening 
coefficient  7,,.  is  often  applied  [9,  11,  12,  13].  The  flow  straightening  coefficient  accounts 
for  the  influence  of  the  hull,  and  typically  has  a  value  of  less  than  one,  meaning  that 
the  cross-flow  velocity  is  reduced  due  to  the  presence  of  the  hull.  Lewandowski  [11] 
gives  the  following  equation  for  estimating  the  flow  straightening  coefficient: 


lr 


1 

1  +  CB 


(40) 


where  Cb  is  the  ship  hull  block  coefficient.  Incorporating  the  effect  of  the  flow  straight¬ 
ening  coefficient,  the  effective  total  velocity  of  the  rudder  relative  to  the  incident  flow 
is  given  by  the  following: 


Vr  = 


(41) 
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where  Fr  is  the  dihedral  angle  of  the  rudder,  as  shown  in  Figure  7.  The  dihedral  angle 
for  a  typical  rudder  is  —90°.  Note  that  terms  with  Tr  in  Equation  (41)  cause  the 
spanwise  component  of  flow  to  be  excluded  when  evaluating  Vr. 


Figure  7:  Dihedral  Angle  I  for  Rudders  and  Other  Foils 

The  flow  angle  of  attack  a  acting  on  the  rudder  includes  contributions  from  both  the 
rudder  deflection  Srudder  and  the  incident  flow.  ShipMo3D  uses  a  rudder  deflection 
convention  of  positive  for  counter-clockwise  rotation  when  viewed  from  inside  the 
ship.  Using  this  convention,  the  total  angle  of  attack  on  a  rudder  is  expressed  as: 


a  =  Sr 


+  OLVr 


The  angle  of  attack  of  the  incident  flow  relative  to  the  x  axis  is  given  by: 


ayr  =  7 r  arctan  —  —  770  )  sin  Tr  —  (  arctan  —  +  775  )  cos  T 


Knowing  the  incident  flow  characteristics,  the  lift  and  drag  forces  acting  on  the  rudder 
can  be  computed  using  Equations  (27)  and  (28),  and  the  forces  acting  on  the  ship 
can  then  be  given  by: 

Fi  =  F[^der  sin  aVr  -  F^r  cos  (44) 

F2  =  {-Ff$der  cos aVr  -  silica)  sin Tr  (45) 

^3  =  {~Fl^r  cos  ayr  -  F™dder  smavr)cosrr  (46) 
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FA  = 

-z  F2  +  y  F3 

(47) 

F5  = 

z  Fi  —  x  F3 

(48) 

F6  = 

-y  Pi  +  x  f2 

(49) 

As  mentioned  previously,  the  above  formulation  gives  forces  in  translating  earth  axes 
assuming  small  amplitude  motions.  Within  ShipMo3D,  large  amplitude  motion  ef¬ 
fects  are  included  by  computing  forces  in  ship-based  axes  and  transforming  values  to 
translating  earth  axes. 

The  above  discussion  of  rudder  forces  does  not  include  lift  forces  due  to  a  seaway.  Due 
to  uncertainties  regarding  interactions  between  the  seaway,  rudder,  and  propeller,  lift 
forces  due  to  the  seaway  are  still  evaluated  without  considering  the  influence  of  the 
propeller,  as  was  originally  done  in  Reference  3.  The  influence  of  this  simplification 
will  likely  be  minor  because  wave-induced  water  particle  velocities  will  typically  be 
small  relative  to  velocities  related  to  ship  maneuvering. 


5  Introduction  of  a 

Proportional-Integral-Derivative  Autopilot 


Due  to  its  widespread  usage  onboard  ships,  a  proportional-integral-derivative  (PID) 
autopilot  has  been  implemented  within  ShipMo3D.  Fossen  [19]  gives  a  comprehensive 
treatment  of  PID  and  other  autopilot  types  for  marine  vehicles. 

The  command  rudder  angle  S}Jdder  given  by  a  PID  autopilot  is  as  follows: 


where  k^  is  the  proportional  gain  for  mode  j,  rjj  is  the  motion  displacement  in  earth- 
fixed  axes  for  mode  j,  rfc-  is  the  command  motion  displacement  for  mode  j,  is 
the  integral  gain  for  mode  j,  T™dder  is  the  integration  duration,  t  is  the  current  time, 
r  is  the  time  delay  for  integration,  kfj  is  the  derivative  gain  for  mode  j,  and  f/:[  is  the 
motion  velocity  in  earth- fixed  axes  for  mode  j.  For  a  given  command  rudder  angle 
^ru,dder ,  w}1jc}1  can  fog  provided  by  the  autopilot  or  by  a  helmsman,  the  rudder  response 
is  modelled  as  follows: 


5r 


+  2  (s  ujrsudder  5r 


+  w? 


=  ^ 


c 


(51) 


where  is  the  nondimensional  damping  response  constant,  and  Us  is  the  rudder 
frequency  response  constant. 
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6  Turning  Circle  Predictions  for  the  Tanker 
Esso  Osaka 


Turning  circle  predictions  have  been  performed  for  the  Esso  Osaka,  an  oil  tanker  that 
has  been  the  subject  of  extensive  maneuvering  studies  in  the  open  literature  [9,  20]. 
Table  1  gives  main  particulars  for  the  Esso  Osaka.  ShipMo3D  predictions  for  the  Esso 
Osaka  were  originally  presented  in  Reference  4. 

Table  1:  Main  Particulars  for  the  Esso  Osaka  During  Maneuvering  Trials 

Length,  L  (between  perpendiculars)  325  m 

Beam,  B  53  m 

Midships  draft,  Tmid  21.73  m 

Trim  by  stern,  tstern  0.0  m 

Displacement,  A  319,400  tonnes 

Longitudinal  CG  forward  of  midship  10.3  m 


6.1  ShipMo3D  Models  of  the  Esso  Osaka 

Two  models  of  the  Esso  Osaka  were  developed  using  the  ShipMo3D  applications 
SM3DPanclHull,  SM3DRadDif,  and  SM3DBuildShip,  which  are  described  in  Refer¬ 
ence  6.  The  first  ShipMo3D  model  of  the  Esso  Osaka  uses  experimental  hull  maneu¬ 
vering  force  coefficients,  while  the  second  model  uses  predicted  hull  maneuvering  force 
coefficients  based  on  Inoue  et  al.  [10].  Both  of  the  ship  models  include  modifications 
to  evaluation  of  hull  maneuvering  forces  and  rudder  forces  described  earlier  in  this 
report.  Figure  8  shows  a  ShipMo3D  model  of  the  Esso  Osaka. 


The  ship  model  based  on  experimental  hull  maneuvering  coefficients  uses  values  for 
a  bare  hull  model  reported  in  Table  5.2  of  the  ITTC  Esso  Osaka  Committee  report 
[9].  Some  of  the  nonlinear  hull  force  coefficients  in  the  ITTC  report  have  a  different 
format  to  those  developed  by  Inoue  et  al.  that  are  used  within  ShipMo3D.  Where 
differences  occur,  equivalent  ShipMo3D  coefficients  can  be  obtained  as  follows: 
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(52) 

(53) 

(54) 

(55) 


where  Y„'„, 
ITTC,  t>' 


Y']rr ,  Y'rrr ,  and  N'r^  are  nonlinear  hull  force  coefficients  reported  by  the 
is  the  nominal  non-dimensional  sway  velocity  for  matching  nonlinear 
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Figure  8:  ShipMo3D  Model  of  Esso  Osaka  Including  Panelled  Wet  Hull,  Propeller, 
and  Rudder 

forces  from  the  fnoue  et  al.  and  fTTC  formulations,  and  r'nom  is  the  nominal  non- 
dimensional  yaw  velocity  for  matching  nonlinear  forces.  The  ShipMo3D  coefficients 
were  obtained  using  turning  circle  characteristics  for  the  Esso  Osaka  with  an  initial 
speed  of  10  knots,  resulting  in  an  equilibrium  condition  with  a  final  total  speed  of  3.5 
knots,  yaw  rate  of  0.3  degrees  per  second,  and  drift  angle  of  15°,  giving  \v'nom\  =  0.262 

and  KoJ  =  °-944- 

For  the  linear  sway-yaw  force  coefficient,  Table  5.2  of  Reference  9  gives  Y'r  —  m'x , 
which  includes  the  nondimensional  surge  added  mass.  Radiation  computations  with 
ShipMo3D  indicate  that  the  nondimensional  surge  added  mass  for  the  Esso  Osaka 
has  a  value  of  0.016  at  zero  forward  speed  and  zero  motion  frequency,  which  has  been 
considered  when  evaluating  an  input  Yx  value  for  ShipMo3D. 

When  using  experimental  hull  coefficients  as  input  to  ShipMo3D,  they  should  be 
referenced  to  the  ship  centre  of  gravity  (CG).  Table  5.2  of  Reference  9  gives  hull 

coefficients  based  on  motions  at  the  CG  and  forces  relative  to  midships.  To  obtain 

forces  relative  to  the  CG  for  ShipMo3D,  the  following  corrections  have  been  applied: 

N'V(CG,  CG)  =  N'v(mid,CG)  +  ^Y'(mid,CG)  (56) 

N’t(CG,CG)  =  N'r(mul,CG)  +  ^  Y^(mid,CG)  (57) 

(58) 

where  N'V(CG,CG )  and  N'r(CG ,CG)  are  hull  force  coefficients  with  both  forces 
and  motions  referenced  to  the  CG,  N'v{mid,CG ),  Y^(mid,CG),  N'r{mid,CG ),  and 
Y'r{mid,CG )  are  hull  force  coefficients  with  forces  referenced  to  midships  and  mo¬ 
tions  referenced  to  the  CG,  and  xmid  is  the  x  coordinate  of  midships  relative  to  the 
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CG.  Note  that  the  sway  force  terms  Y^(mid,CG)  and  Y^(mid,CG)  require  no  cor¬ 
rections  to  obtain  values  with  forces  referenced  to  the  ship  CG.  For  nonlinear  hull 
force  terms,  it  is  assumed  that  corrections  of  force  coefficients  to  be  relative  to  the 
CG  are  negligible. 

Table  2  shows  generally  good  agreement  between  predicted  and  experimental  hull 
force  coefficients.  The  experimental  coefficients  in  Table  2  are  given  for  both  forces 
and  motions  relative  to  the  CG.  For  the  predicted  hull  force  coefficients  based  on 
Inoue  et  al.,  no  corrections  have  been  applied  to  transfer  force  and  motion  locations 
from  midships  to  the  CG  because  the  magnitude  of  such  corrections  would  be  small 
relative  to  the  uncertainties  in  predicted  hull  force  coefficients.  When  comparing  the 
hull  force  coefficients  in  Table  2,  the  most  notable  differences  are  for  the  terms  Y'r ^ 
and  N'vr 2,  which  both  have  predicted  values  of  0.0. 

Table  2:  Esso  Osaka  Hull  Force  Coefficients  from  Experiments  [9]  and  Predicted  Using 
Inoue  et  al.  [10] 


Experiments 

Predicted 

Sway-sway 

Y' 

V 

-0.383 

-0.400 

y> 

1v\v\ 

-0.276 

-0.360 

Y> 

1v\r\ 

-0.242 

-0.317 

Sway-yaw 

Y' 

r 

0.098 

0.105 

Y' 

1  r\r\ 

-0.011 

0.0 

Yaw-sway 

K 

-0.135 

-0.134 

Kr> 

0.024 

0.000 

Yaw-yaw 

K 

-0.051 

-0.054 

Ky 

-0.297 

-0.200 

K\rl 

-0.017 

-0.012 

The  resistance  and  propulsion  characteristics  for  the  Esso  Osaka  were  assumed  to 
be  the  same  as  from  the  previous  ShipMo3D  maneuvering  study  [4],  Table  3  gives 
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Table  3:  Estimated  Resistance  Coefficients  for  Esso  Osaka 


Speed 

(knots) 

2.0 

16.0 

18.0 

19.6 

21.2 

22.9 


Resistance  coefficient 
CDx  =  F  /  (l  /  2  pU2Aw) 
0.00306 
0.00306 
0.00308 
0.00328 
0.00370 
0.00426 


estimated  resistance  coefficients  based  on  the  Series  60  hull  with  a  block  coefficient  of 
0.80  [21],  which  has  a  geometry  very  similar  to  the  Esso  Osaka.  Based  on  data  given 
in  Reference  22,  the  wake  fraction  wprop  was  estimated  to  be  0.352  and  the  thrust 
deduction  coefficient  tprop  was  estimated  to  be  0.2.  The  Esso  Osaka  has  a  propeller 
diameter  of  9.1  m.  The  following  thrust  coefficient  curve  was  obtained  by  scaling  a 
reference  curve  to  match  the  reported  full-scale  speed  of  10  knots  with  a  propeller 
speed  of  51  RPM. 


Kr 


0.394  -  0.197  Jprop 


0.148  Jprop 


(59) 


Table  4  gives  rudder  dimensions,  autopilot  parameters,  and  hydrodynamic  coeffi¬ 
cients.  The  rudder  wake  fraction  and  flow  straightening  coefficient  are  based  on  re¬ 
ported  values  in  Reference  9.  The  rudder-propeller  interaction  coefficient  is  based  on 
the  ratio  of  the  propeller  diameter  to  the  rudder  span,  and  represents  an  estimate  of 
the  fraction  of  the  rudder  within  the  propeller  slipstream. 

6.2  Comparisons  of  Predicted  and  Measured  Turning 
Circles 

Figures  9  to  14  show  comparisons  of  predicted  and  measured  turning  circles  for  the 
Esso  Osaka.  When  considering  the  measured  turning  circles,  it  should  be  noted  that 
scatter  is  greater  for  the  lower  speed  of  7.7  knots,  likely  due  to  the  greater  susceptibil¬ 
ity  to  environment  forces  at  lower  ship  speed.  During  the  initial  stages  of  each  turning 
circle,  both  sets  of  predictions  agree  closely  with  the  measured  results.  Differences 
among  the  measured  results  and  two  sets  of  predictions  increase  as  the  turn  radius 
decreases,  which  is  expected  because  nonlinear  forces,  which  are  difficult  to  predict, 
become  more  important.  Of  the  two  sets  of  predictions,  the  predictions  based  on  the 
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Table  4:  Rudder  Dimensions,  Autopilot  Parameters,  and  Hydrodynamic  Coefficients 
for  Esso  Osaka 


Rudder  span  srudder 

13.85  m 

Rudder  chord  length  crudder 

9.00  m 

Maximum  rudder  deflection  S^dder 

35  degrees 

Maximum  rudder  velocity  S7n^er 

3  deg/s 

Response  natural  frequency  c o$ 

3  rad/s 

Response  damping 

0.85 

Rudder  wake  fraction  wr 

0.0 

Rudder  flow  straightening  coefficient  yr 

0.4 

Rudder-propeller  interaction  coefficient  crudder~prop 

0.66 

experimental  hull  force  coefficients  appear  to  be  better  than  the  predictions  based  on 
hull  coefficients  estimated  based  on  Inoue  et  ah,  as  would  be  expected.  The  estimated 
hull  coefficients  give  lower  speed  and  higher  yaw  rates,  resulting  in  noticeably  smaller 
final  turning  radius. 
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xf  (m) 


□  Full-scale  trials 


Predictions,  experimental  hull  coefficients 


Figure  9:  Esso  Osaka  Trajectory  During  Turning  Circle,  Initial  Speed  of  10  knots 
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□  Full-scale  trials 


Predictions,  experimental  hull  coefficients 


Figure  10:  Esso  Osaka  Total  Speed  During  Turning  Circle,  Initial  Speed  of  10  knots 


□  Full-scale  trials 


Figure  11:  Esso  Osaka  Yaw  Rate  During  Turning  Circle,  Initial  Speed  of  10  knots 
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□  Full-scale  trials 


-  Predictions,  experimental  hull  coefficients 


Predictions,  Inoue  hull  coefficients 


yf  (m) 


Figure  12:  Esso  Osaka  Trajectory  During  Turning  Circle,  Initial  Speed  of  7. 7  knots 
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xf  (m) 


□  Full-scale  trials 


Figure  13:  Esso  Osaka  Total  Speed  During  Turning  Circle,  Initial  Speed  of  7.7  knots 


t(s) 


Figure  14:  Esso  Osaka  Yaw  Rate  During  Turning  Circle,  Initial  Speed  of  7.7  knots 


DRDC  Atlantic  TM  2008-162 


23 


7  Conclusions 


Improvements  have  been  made  to  evaluation  of  hull  maneuvering  forces  and  rudder 
forces  within  ShipMo3D.  ShipMo3D  also  now  has  a  new  capability  for  modelling  PID 
autopilots.  Turning  circle  predictions  for  the  tanker  Esso  Osaka  show  good  agreement 
with  sea  trials  data.  Turning  circle  predictions  based  on  experimental  hull  maneu¬ 
vering  coefficients  are  more  accurate  than  predictions  based  on  hull  maneuvering 
coefficients  estimated  using  the  regression  method  of  Inoue  et  al.  [10]. 
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Symbols  and  Abbreviations 


M 

[Abk~i] 

| jsjoil-i 

dA 

^j^rudder 
^ skeg-i 


'  hull  , 
jk  /dx 


B 

[B] 

[5^] 

[Bfoil-ij 

^jghull— manj 
]^phull— rad^ 

^skeg-ij 

ryhull—man 
n22V 
jDhull—man 
n22V\r'\ 
hull— man 


'22V\v'\ 
lhull—man 
>26V 

^ hull— man 
>2QV\r'\ 
ryhull—man 
n62V 
r>hull—man 
n62V\r'2\ 

^ hull— man 

>mv 

^ hull— man 
;66y|r'| 

TDhull—man 

n6QV\v'2\ 

-ihull— rad— stab 
>22  U 

-ihull— rad— stab 
>2  QU 
r>hull— rad— stab 
n62U 

TDhull— rad— stab 

nmu 

CB 

Cdx 

ibk—il 


B'2 

bL 


Bq 

Bl 


Bo, 

Bl 


[Cb 


added  mass  matrix 
added  mass  matrix  of  bilge  keel  i 
added  mass  matrix  of  foil  i 
hull  added  mass  matrix 

term  based  on  x  derivatives  of  added  mass  potentials 

added  mass  matrix  of  rudder  i 

added  mass  matrix  of  skeg  i 

effective  aspect  ratio  of  foil 

ship  beam 

damping  matrix 

damping  matrix  of  bilge  keel  i 

damping  matrix  of  foil  i 

hull  maneuvering  damping  matrix 

frequency- independent  damping  matrix  due  to  potential  flow 

hull  viscous  damping  matrix 

damping  matrix  of  skeg  i 

sway-sway  linear  maneuvering  coefficient 

sway-sway  maneuvering  coefficient  for  yaw  nonlinearity 

sway-sway  maneuvering  coefficient  for  sway  nonlinearity 

sway-yaw  linear  maneuvering  coefficient 

sway-sway  maneuvering  coefficient  for  yaw  nonlinearity 

yaw-sway  linear  maneuvering  coefficient 

yaw-sway  maneuvering  coefficient  for  yaw  nonlinearity 

yaw-yaw  linear  maneuvering  coefficient 

yaw-yaw  maneuvering  coefficient  for  yaw  nonlinearity 

yaw-yaw  maneuvering  coefficient  for  sway  nonlinearity 

sway-sway  radiation  damping  term 

sway-yaw  radiation  damping  term 

yaw-sway  radiation  damping  term 

yaw-yaw  radiation  damping  term 

ship  block  coefficient 

hull  resistance  coefficient 

stiffness  force  matrix  for  bilge  keel  i 
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[(Jbuoyj 
(jdrag 
| (jfoil-i^ 

(jr  udder  —prop 


3 


]^jhull— manj 
(Jlift 
]^jskeg— zj 

cth 

CG 

c 

Dprop 

Fj 

{FD} 

{F1} 

|  pi-bk-iy 

|  pl-foil-iy 

r pl—hull\ 

|  pi— skeg— i\ 

T^drag 
rlift 
rprudder 
r  lift 

jpprop 

{FruJJer-,} 

Fi 

Jprop 

kt 

^phullj 

hD 

KSj 

k1 

KSj 

kp 

KSj 

L 

[M] 

m'x 

NL 


hull  buoyancy  stiffness  matrix 
foil  drag  force  coefficient 
stiffness  force  matrix  for  foil  i 

rudder-propeller  interaction  coefficient  from  propeller  j 

hull  maneuvering  force  stiffness  matrix 

foil  lift  force  coefficient 

stiffness  force  matrix  for  skeg  % 

propeller  thrust  loading  coefficient 

centre  of  gravity 

mean  chord  length 

hull  frequency-independent  stiffness  due  to  potential  flow 
propeller  diameter 

force  component  for  mode  j  in  translating  earth  axes 

wave  diffraction  force  vector 

incident  wave  excitation  force  vector 

incident  wave  force  vector  on  bilge  keel  i 

incident  wave  force  vector  on  foil  i 

incident  wave  force  vector  on  hull 

incident  wave  force  vector  on  skeg  i 

rudder  drag  force 

rudder  lift  force 

propeller  thrust 

non-inertial  force  vector  on  rudder  i 

force  component  for  mode  j  in  ship-based  axes 

propeller  advance  ratio 

propeller  thrust  coefficient 

retardation  function  matrix 

autopilot  derivative  gain  for  mode  j 

autopilot  integral  gain  for  mode  j 

autopilot  proportional  gain  for  mode  j 

ship  length  between  perpendiculars 

ship  inertia  matrix 

nondimensional  surge  added  mass 

non-dimensional  yaw-yaw  maneuvering  force  coefficient 
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N' 

V 


^ slip 

rOO 

'  slip 

s 


s 

T1 

mid 

t 

t  stern 

u 


ur ,  vr,  wr 

V'slip—j 


u 


pot 

slip 


U 


slip 


v 

Vr 


w 


prop 

xf,Vf 

xf,yf 

Y’r 

Y’ 


Y' 

vrr 


Y' 

^  vvv 


%lift 

a 

Y 

^ rudder 


non-dimensional  yaw-yaw  maneuvering  force  coefficient  with  yaw  non¬ 
linearity 

non-dimensional  yaw-sway  maneuvering  force  coefficient 
non-dimensional  yaw  velocity 
nominal  non-dimensional  yaw  velocity 
slipstream  radius 

slipstream  radius  far  aft  of  propeller 
foil  area 
foil  span 

draft  at  midships 
time 

trim  by  stern 

ship  forward  speed 

water  relative  velocity  components 

slipstream  velocity  aft  of  propeller  j 

slipstream  velocity  aft  of  propeller  based  on  potential  flow 

flow  velocity  far  aft  of  propeller 

total  ship  speed  in  horizontal  plane 

incident  flow  velocity  on  rudder 

non-dimensional  sway  velocity 

nominal  non-dimensional  sway  velocity 

propeller  wake  fraction 

horizontal  plane  coordinates  in  earth-fixed  axes 
horizontal  plane  accelerations 

non-dimensional  sway-yaw  maneuvering  force  coefficient 

non-dimensional  sway-yaw  maneuvering  force  coefficient  with  yaw  non¬ 
linearity 

non-dimensional  sway-sway  maneuvering  force  coefficient 
non-dimensional  sway-sway  maneuvering  force  coefficient  with  yaw  non¬ 
linearity 

non-dimensional  sway-sway  maneuvering  force  coefficient  with  sway 
nonlinearity 

longitudinal  location  of  lift  force 
flow  angle  of  attack 
rudder  flow  straightening  constant 
rudder  deflection  angle 
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$maxer  maximum  rudder  deflection  angle 
firudder  rudder  deflection  velocity 
S™ax£r  maximum  rudder  deflection  velocity 
fir-udder-  ru(](]er  deflection  acceleration 
firudder  command  rudder  angle 
Q  rudder  control  damping  coefficient 

Tjj  motion  displacement  in  translating  earth  axes  for  mode  j 

ijj  acceleration  in  translating  earth  axes  for  mode  j 

r)'-  motion  displacement  in  earth-fixed  axes  for  mode  j 

rfC]  command  motion  displacement  in  earth-fixed  axes  for  mode  j 

ijj  motion  velocity  in  earth-fixed  axes  for  mode  j 

p  water  density 

Tmaxer  autopilot  integration  duration 

X  yaw  velocity  in  earth-fixed  axes 

X  yaw  acceleration  in  earth-fixed  axes 

uje  encounter  frequency 

0J5  rudder  control  natural  frequency 

A  ship  mass  displacement 
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